The L1157-B1 astrochemical laboratory: testing the origin of DCN

L1157-B1 is the brightest shocked region of the large-scale molecular outflow, considered the prototype of chemically rich outflows, being the ideal laboratory to study how shocks affect the molecular gas. Several deuterated molecules have been previously detected with the IRAM 30m, most of them formed on grain mantles and then released into the gas phase due to the shock. We aim to observationally investigate the role of the different chemical processes at work that lead to formation the of DCN and test the predictions of the chemical models for its formation. We performed high-angular resolution observations with NOEMA of the DCN(2-1) and H13CN(2-1) lines to compute the deuterated fraction, Dfrac(HCN). We detected emission of DCN(2-1) and H13CN(2-1) arising from L1157-B1 shock. Dfrac(HCN) is ~4x10$^{-3}$ and given the uncertainties, we did not find significant variations across the bow-shock. Contrary to HDCO, whose emission delineates the region of impact between the jet and the ambient material, DCN is more widespread and not limited to the impact region. This is consistent with the idea that gas-phase chemistry is playing a major role in the deuteration of HCN in the head of the bow-shock, where HDCO is undetected as it is a product of grain-surface chemistry. The spectra of DCN and H13CN match the spectral signature of the outflow cavity walls, suggesting that their emission result from shocked gas. The analysis of the time dependent gas-grain chemical model UCL-CHEM coupled with a C-type shock model shows that the observed Dfrac(HCN) is reached during the post-shock phase, matching the dynamical timescale of the shock. Our results indicate that the presence of DCN in L1157-B1 is a combination of gas-phase chemistry that produces the widespread DCN emission, dominating in the head of the bow-shock, and sputtering from grain mantles toward the jet impact region.Comment: Accepted for publication in A&A. 7 pages, 5 Figures, 1 Tabl

Temperature and kinematics of protoclusters with intermediate and high-mass stars: the case of IRAS 05345+3157

We have mapped at small spatial scales the temperature and the velocity field in the protocluster associated with IRAS 05345+3157, which contains both intermediate-/high-mass protostellar candidates and starless condensations, and is thus an excellent location to investigate the role of massive protostars on protocluster evolution. We observed the ammonia (1,1) and (2,2) inversion transitions with the VLA. Ammonia is the best thermometer for dense and cold gas, and the observed transitions have critical densities able to trace the kinematics of the intracluster gaseous medium. The ammonia emission is extended and distributed in two filamentary structures. The starless condensations are colder than the star-forming cores, but the gas temperature across the whole protocluster is higher (by a factor of ~1.3-1.5) than that measured typically in both infrared dark clouds and low-mass protoclusters. The non-thermal contribution to the observed line broadening is at least a factor of 2 larger than the expected thermal broadening even in starless condensations, contrary to the close-to-thermal line widths measured in low-mass quiescent dense cores. The NH3-to-N2H+ abundance ratio is greatly enhanced (a factor of 10) in the pre--stellar core candidates, probably due to freeze-out of most molecular species heavier than He. The more massive and evolved objects likely play a dominant role in the physical properties and kinematics of the protocluster. The high level of turbulence and the fact that the measured core masses are larger than the expected thermal Jeans masses indicate that turbulence likely was an important factor in the initial fragmentation of the parental clump.Comment: 13 pages (with Appendix), 11 figure

Thermal Jeans fragmentation within 1000 AU in OMC-1S

We present subarcsecond 1.3 mm continuum ALMA observations towards the Orion Molecular Cloud 1 South (OMC-1S) region, down to a spatial resolution of 74 AU, which reveal a total of 31 continuum sources. We also present subarcsecond 7 mm continuum VLA observations of the same region, which allow to further study fragmentation down to a spatial resolution of 40 AU. By applying a Mean Surface Density of Companions method we find a characteristic spatial scale at ~560 AU, and we use this spatial scale to define the boundary of 19 `cores' in OMC-1S as groupings of millimeter sources. We find an additional characteristic spatial scale at ~2900 AU, which is the typical scale of the filaments in OMC-1S, suggesting a two-level fragmentation process. We measured the fragmentation level within each core and find a higher fragmentation towards the southern filament. In addition, the cores of the southern filament are also the densest (within 1100 AU) cores in OMC-1S. This is fully consistent with previous studies of fragmentation at spatial scales one order of magnitude larger, and suggests that fragmentation down to 40 AU seems to be governed by thermal Jeans processes in OMC-1S.Comment: Accepted to Ap

The B1 shock in the L1157 outflow as seen at high spatial resolution

We present high spatial resolution (750 AU at 250 pc) maps of the B1 shock in the blue lobe of the L1157 outflow in four lines: CS (3-2), CH3OH (3_K-2_K), HC3N (16-15) and p-H2CO (2_02-3_01). The combined analysis of the morphology and spectral profiles has shown that the highest velocity gas is confined in a few compact (~ 5 arcsec) bullets while the lowest velocity gas traces the wall of the gas cavity excavated by the shock expansion. A large velocity gradient model applied to the CS (3-2) and (2-1) lines provides an upper limit of 10^6 cm^-3 to the averaged gas density in B1 and a range of 5x10^3< n(H2)< 5x10^5 cm^-3 for the density of the high velocity bullets. The origin of the bullets is still uncertain: they could be the result of local instabilities produced by the interaction of the jet with the ambient medium or could be clump already present in the ambient medium that are excited and accelerated by the expanding outflow. The column densities of the observed species can be reproduced qualitatively by the presence in B1 of a C-type shock and only models where the gas reaches temperatures of at least 4000 K can reproduce the observed HC3N column density.Comment: 13 pages, 12 figure

Deuteration as an evolutionary tracer in massive-star formation

Theory predicts, and observations confirm, that the column density ratio of a molecule containing D to its counterpart containing H can be used as an evolutionary tracer in the low-mass star formation process. Since it remains unclear if the high-mass star formation process is a scaled-up version of the low-mass one, we investigated whether the relation between deuteration and evolution can be applied to the high-mass regime. With the IRAM-30m telescope, we observed rotational transitions of N2D+ and N2H+ and derived the deuterated fraction in 27 cores within massive star-forming regions understood to represent different evolutionary stages of the massive-star formation process. Results. Our results clearly indicate that the abundance of N2D+ is higher at the pre-stellar/cluster stage, then drops during the formation of the protostellar object(s) as in the low-mass regime, remaining relatively constant during the ultra-compact HII region phase. The objects with the highest fractional abundance of N2D+ are starless cores with properties very similar to typical pre-stellar cores of lower mass. The abundance of N2D+ is lower in objects with higher gas temperatures as in the low-mass case but does not seem to depend on gas turbulence. Our results indicate that the N2D+-to-N2H+ column density ratio can be used as an evolutionary indicator in both low- and high-mass star formation, and that the physical conditions influencing the abundance of deuterated species likely evolve similarly during the processes that lead to the formation of both low- and high-mass stars.Comment: Accepted by A&AL, 4 pages, 2 figures, 2 appendices (one for Tables, one for additional figures

Broad N2H+ emission towards the protostellar shock L1157-B1

We present the first detection of N2H+ towards a low-mass protostellar outflow, namely the L1157-B1 shock, at about 0.1 pc from the protostellar cocoon. The detection was obtained with the IRAM 30-m antenna. We observed emission at 93 GHz due to the J = 1-0 hyperfine lines. The analysis of the emission coupled with the HIFI CHESS multiline CO observations leads to the conclusion that the observed N2H+(1-0) line originates from the dense (> 10^5 cm-3) gas associated with the large (20-25 arcsec) cavities opened by the protostellar wind. We find a N2H+ column density of few 10^12 cm-2 corresponding to an abundance of (2-8) 10^-9. The N2H+ abundance can be matched by a model of quiescent gas evolved for more than 10^4 yr, i.e. for more than the shock kinematical age (about 2000 yr). Modelling of C-shocks confirms that the abundance of N2H+ is not increased by the passage of the shock. In summary, N2H+ is a fossil record of the pre-shock gas, formed when the density of the gas was around 10^4 cm-3, and then further compressed and accelerated by the shock.Comment: ApJ, in pres

Dense gas in IRAS 20343+4129: an ultracompact HII region caught in the act of creating a cavity

The intermediate- to high-mass star-forming region IRAS 20343+4129 is an excellent laboratory to study the influence of high- and intermediate-mass young stellar objects on nearby starless dense cores, and investigate for possible implications in the clustered star formation process. We present 3 mm observations of continuum and rotational transitions of several molecular species (C2H, c-C3H2, N2H+, NH2D) obtained with the Combined Array for Research in Millimetre-wave Astronomy, as well as 1.3 cm continuum and NH3 observations carried out with the Very Large Array, to reveal the properties of the dense gas. We confirm undoubtedly previous claims of an expanding cavity created by an ultracompact HII region associated with a young B2 zero-age main sequence (ZAMS) star. The dense gas surrounding the cavity is distributed in a filament that seems squeezed in between the cavity and a collimated outflow associated with an intermediate-mass protostar. We have identified 5 millimeter continuum condensations in the filament. All of them show column densities consistent with potentially being the birthplace of intermediate- to high-mass objects. These cores appear different from those observed in low-mass clustered environments in sereval observational aspects (kinematics, temperature, chemical gradients), indicating a strong influence of the most massive and evolved members of the protocluster. We suggest a possible scenario in which the B2 ZAMS star driving the cavity has compressed the surrounding gas, perturbed its properties and induced the star formation in its immediate surroundings.Comment: 17 pages, 13 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society (Main Journal

Mid-J CO Emission in Nearby Seyfert Galaxies

We study for the first time the complete sub-millimeter spectra (450 GHz to 1550 GHz) of a sample of nearby active galaxies observed with the SPIRE Fourier Transform Spectrometer (SPIRE/FTS) onboard Herschel. The CO ladder (from Jup = 4 to 12) is the most prominent spectral feature in this range. These CO lines probe warm molecular gas that can be heated by ultraviolet photons, shocks, or X-rays originated in the active galactic nucleus or in young star-forming regions. In these proceedings we investigate the physical origin of the CO emission using the averaged CO spectral line energy distribution (SLED) of six Seyfert galaxies. We use a radiative transfer model assuming an isothermal homogeneous medium to estimate the molecular gas conditions. We also compare this CO SLED with the predictions of photon and X-ray dominated region (PDR and XDR) models.Comment: Proceedings of the Torus Workshop 2012 held at the University of Texas at San Antonio, 5-7 December 2012. C. Packham, R. Mason, and A. Alonso-Herrero (eds.); 6 pages, 3 figure

Iron and Nickel spectral opacity calculations in conditions relevant for pulsating stellar envelopes and experiments

Seismology of stars is strongly developing. To address this question we have formed an international collaboration OPAC to perform specific experimental measurements, compare opacity calculations and improve the opacity calculations in the stellar codes [1]. We consider the following opacity codes: SCO, CASSANDRA, STA, OPAS, LEDCOP, OP, SCO-RCG. Their comparison has shown large differences for Fe and Ni in equivalent conditions of envelopes of type II supernova precursors, temperatures between 15 and 40 eV and densities of a few mg/cm3 [2, 3, 4]. LEDCOP, OPAS, SCO-RCG structure codes and STA give similar results and differ from OP ones for the lower temperatures and for spectral interval values [3]. In this work we discuss the role of Configuration Interaction (CI) and the influence of the number of used configurations. We present and include in the opacity code comparisons new HULLAC-v9 calculations [5, 6] that include full CI. To illustrate the importance of this effect we compare different CI approximations (modes) available in HULLAC-v9 [7]. These results are compared to previous predictions and to experimental data. Differences with OP results are discussed.Comment: 4 pages, 3 figures, conference Inertial Fusion Sciences and Applications, Bordeaux, 12th to 16th September 2011; EPJ web of Conferences 201